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Guided  ion  beam  kinetic  energy  thresholds  in  the  ion- 
molecule  reactions  M++H2->MH++H  where  M+  is  a  closed-shell 

atomic  ion  B+,  Al+,  or  Ga+,  were  found  to  exceed  by  0.4  to 
c.a.  5  eV  the  thermodynamic  energy  requirements  (or  the 
theoretically  computed  barrier  heights)  for  these  reactions. 

In  addition,  the  formation  of  MD+  occurs  at  a  significantly 
lower  threshold  than  MH+  when  M+  reacts  with  HD.  Moreover, 
the  measured  reaction  cross-sections  for  the  production  of 
MH+  or  MD+  product  ions  are  very  small  (10-17  to  10~20 
cm2), being  largest  for  B+  and  smallest  for  Ga+.  A  previous 
paper  from  this  group  proposed  that  collisional-to-internal 
energy  transfer  is  the  rate-limiting  step  for  this  class  of 
reactions.  It  also  suggested,  based  on  a  dynamical  resonance 
picture,  that  collisions  occurring  at  or  near  C2v  symmetry  are 
more  effective  than  other  collisions  even  though  C2v 
geometries  provide  no  lower  potential  energy  barriers  than 
others.  By  examining  the  collision  paths  characteristic  of 
flux  early  in  the  bimolecular  collision  and  searching  for 
geometries  along  such  paths  where  collision- to- internal 
energy  transfer  is  optimal,  our  earlier  efforts  predicted 
reaction  thresholds  in  reasonable  agreement  with  the 
(previously  perplexing)  experimental  data.  In  the  present 
work,  we  introduce  a  model  Hamiltonian  whose  classical  and 
quantum  dynamics  we  apply  to  the  M+  +  H2,  D2,  HD  reactive 
collisions.  We  calculate  the  classical  collisional-to- 
internal  energy  transfer  cross-section  and  find  energy 
transfer  thresholds  that  resemble  the  experimental  reaction 
thresholds  but  whose  isotopic  mass  trends  are  not  entirely 
consistent  with  experiment.  We  then  use  a  Green  Function 
method  and  a  local  quadratic  approximation  to  the  potential 
surface  to  obtain  analytical  expressions  for  the  isotopic 
mass  dependences  of  the  collisional-to-vibrational  energy 
transfer  and  for  the  subsequent  fragmentation  of  the  three- 
atom  system.  Finally,  we  analyze  the  origin  of  the  threshold 
energy  asymmetry  in  the  M++HD  reactions. 


X.  Introduction 


Guided  ion  beam  measurements  of  the  cross-sections1  for 

the  production  of  MH+  and  MD*  product  ions  in  reactions  of 

closed-shell  ^-S  B+,  Al+,  and  Ga+  ions  (denoted  M*)  with 

closed-shell  H2,  D2  and  HD  displayed  features  that 

required  further  interpretation: 

•  The  apparent  thresholds  (i.e.  the  collision  kinetic 
energies  where  product  MH+  or  MD+  ions  are  first  formed) 
exceed  the  minimum  thermodynamic  energy  requirements  by 
significant  amounts  (e.g.,  by  up  to  5  eV  for  Ga+)  .  Two 
examples  of  this  data  are  shown  in  Fig.  1. 

•  In  experiments  with  HD,  MD+  formation  displays  a  lower 
energy  threshold  than  MH+  (see  Fig.  1)  . 

•  The  cross  sections  are  small  (lO-^-7  to  10~20  cm2),  and  are 
smallest  for  Ga+  and  largest  for  B+. 


Impulsive,  statistical,  and  spectator-stripper  models2 
do  not  succeed  in  rationalizing  the  unexpectedly  high 
threshold  energies  or  the  magnitudes  of  the  MD+/MH+  threshold 
energy  asymmetries.  For  example,  an  impulsive  model 

decomposes  the  collision  kinetic  energy  T  for  the  M++HD  case 

1  2  ,  . 
into  components  =  —  T  and  TMD  =  T  that  give  the  kinetic 

energy  of  M+  relative  to  the  H  and  D  atoms,  respectively. 

Such  a  model  then  predicts  that  MD+  formation  can  occur  at  a 


1 


lower  total  collision  energy  T  because  2/3  of  this  energy  is 
available  to  the  M-D  coordinate.  In  particular,  the  model 


predicts  that  the  threshold  kinetic  energy  for  MD+  formation 
^  *  that  for  MH+  formation  in  the  M*+HD 


should  be 


u1' 
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case;  this  quantitative  prediction  is  not  seen  in  the 
experimental  data.  In  contrast,  a  previous  work3  by  one  of 
the  authors  and  collaborators  proposed  that  entrance-channel 
transfer  of  M+-H2  relative  translational  energy  to  H2 
vibrational  energy  is  the  rate-limiting  step  in  this  class  of 
reactions . 


By  considering  collision  paths  characteristic  of  early 
reactant  flux,  in  which  the  H2  D2  or  HD  internuclear  distance 
is  essentially  undisturbed  from  its  equilibrium  value 
(because  the  experiments  involve  room  temperature  Hydrogen 
gas)  and  searching,  along  these  paths,  for  geometries  at 
which  energy  transfer  is  predicted  by  a  resonance  condition 
to  be  optimal,  ref.  3  made  predictions  of  reaction  thresholds 
reasonably  in  agreement  with  the  experimental  findings.  It 
is  the  purpose  of  the  present  work  to  make  more  quantitative 
the  energy- transfer  rate-determining-step  picture  and  to 
carry  out  classical  simulations  and  quantal  dynamical 
analyses  on  a  model  potential  surface  fit  to  the  ab  initio 
data  of  ref.  3  to  gain  further  insight  into  this  class  of 
reactions. 

We  approach  the  problem  in  the  following  manner: 
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1.  We  fit  the  same  ab-initio  C2v  and  near  C2v  potential 
energy  surfaces  employed  in  ref.  3  (see  Fig.  2)  to 
extract  strength  and  range  parameters  characteristic  of 
the  repulsive  portion  of  the  M+-H2  surface  implicated  in 
ref.  3  as  the  regions  near  which  dynamical  resonances  can 
occur.  We  use  the  undisturbed  H-H  potential  to 
characterize  the  H2,  D2,  or  HD  vibration  (because  ref.  3 
showed  the  dynamical  resonances  to  occur  at  geometries 
where  one  of  the  local  normal  mode  frequencies  was  nearly 
equal  to  that  of  the  isolated  hydrogenic  molecule) .  And 
we  use  limited  locally  computed  ab  initio  estimates  of 
the  coupling  strength  between  the  H-H  and  M+-H2 
coordinates . 

2.  We  employ  the  potential  created  in  step  1  to  introduce  a 
two-dimensional  model  Hamiltonian  to  describe  the 
collisional  energy  transfer  entrance -channel  part  of  the 
reaction. 

3.  We  use  an  adaptive-step-size  fifth-order  Runge-Kutta 
method4  to  propagate  classical  trajectories  and  obtain  the 
cross  section  for  the  T-V  energy  transfer.  Trends  in  the 
T-v  energy  transfer  onsets  correlate  with  the 
experimental  thresholds,  and  v=0-»v=l  excitation  is  seen 

to  dominate  near  the  threshold.  However,  the  isotopic 
mass  effects  displayed  in  the  classical  trajectory  data 
are  not  entirely  consistent  witr  the  experimental  results 
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(which  suggests  that  there  is  more  to  the  reaction  cross- 
section  than  T-V  energy  transfer) . 

4.  We  also  analyze  the  quantum  transition  probabilities  for 
the  T-V  process  using  a  Green  Function  method,  to  further 
quantify  the  origin  of  the  high  energy  threshold  and  to 
achieve  analytical  expressions  for  the  threshold  energies 
that  more  clearly  display  the  isotopic  masses  and 
potential  surface  parameters. 

5.  Finally,  we  analyze  energy  transfer  for  the  HD  case,  and 
offer  an  explanation  for  the  origin  of  the  asymmetry  in 
the  MD+  and  MH+  thresholds. 

In  Section.  II  we  develop  the  model  Hamiltonian  and 
describe  the  results  of  classical  trajectory  simulations 
using  it.  In  Section  III,  we  describe  the  quantum  propagator 
approach  and  make  connection  with  the  mode  resonance  picture 
of  ref.  3.  We  also  use  classical  coupled-oscillator  concepts 
to  address  the  threshold  energy  asymmetry  in  the  M++HD  cases. 
Section  IV  contains  an  overview  of  our  findings. 


II.  The  Model  Hamiltonian  and  Classical  Trajectories 

II. a.  The  Model  Potential 

Based  on  the  evidence  detailed  in  ref.  3,  we  assume  that 
near  C2V  collision  geometries,  with  the  H2  moiety  near  its 
equilibrium  internuclear  distance,  are  optimally  effective  at 
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allowing  T-V  energy  transfer  to  occur.  Of  course,  collisions 
occur  at  many  orientations  and  many  impact  parameters  that 
must  be  averaged  over  to  compute  the  cross-section.  However 
we  expect  from  the  earlier  work  that  near-C2v  geometries, 
which  necessarily  have  small  impact  parameters,  will  dominate 
the  T-V  process  because  only  near  such  geometries  does 
dynamical  resonance  occur  between  the  M+  -  H2  collisional  and 
H-H  vibrational  modes.  We  therefore  introduce  two  coordinates 
x  and  y  to  describe  the  BC  +  A  collision  as  shown  in  Fig.  3. 


The  M++H-H  internal  potential  energy  U  is  thus  assumed 
to  be  a  function  of  these  two  coordinates  and  the  following 
two  dimensional  model  Hamiltonian  is  introduced: 


h  d 


h  3 2 


H  - - —  ~r~r - -r-r  +  U(x,  y ) 

2  m  dx.  2  H  <?y‘ 


As  explained  in  the  Introduction,  we  choose  to  express  U(x,y) 

in  terms  of  a  M+-H2  repulsion  which  we  describe  as  v  exp(-ax) , 

1  /  \2 

the  unperturbed  H-H  potential  -  kBC[y  -  yeq )  ,  plus  a  coupling 

2 

which  we  embody  in  a  b  parameter  computed  as  detailed  below: 
U(x ,  y)  =  v  •  exp  (-ax  -  b(y  -  yeq))  +  j  kBC(y  -  yeq  )2 . 


In  this  Hamiltonian,  y  is  the  distance  between  B  and  C,  and  x 
is  the  distance  between  A  and  the  center  of  mass  of  BC. 

The  potential  parameters  a  and  v  are  obtained  by  fitting 
the  ab-initio  computed  data  of  ref. 3  (plus  additional  data 
obtained  by  us)  at  several  values  of  x  but  with  y=yeq»  to  the 
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v  exp (-ax)  functional  form.  Specifically  data  in  the  narrow 
entrance  channels  of  Fig.  2  (r  »  0.7A  R  >  1.0A)  are  used. 

Examples  of  such  fits  are  shown  in  Fig.  4. 

The  b  parameter  describes  the  coupling  between  the 
collision  (x)  and  vibrational  (y)  coordinates.  Here,  xA  is  a 

representative  value  of  x  in  the  region  where  the  model  of 
ref.  3  suggests  energy  transfer  should  occur  (see  Fig.  2), 
U(xA,y)  is  computed  at  several  values  of  y  (ranging  ca.  ±  0.4 

Bohr  from  yeq)  in  a  fully  ab  initio  manner.  The  quantity  In 
[U(xA,y)  -  j^bc  (y-yeq)2!  is  then  plotted  vs  y.  According  to 

our  model  potential,  this  logarithmic  function  should  reduce 
to  In  v  -  axA  -  b(y-yeq)  ,  so  the  slope  of  such  a  plot  should 
give  b.  Examples  of  such  plots  are  given  in  Fig.  5. 


The  values  of  v,  a  and  b  thus  obtained  and  used  for  the 
reactions  studied  here  are  given  in  Table  1. 


The  appropriate  effective  masses  appearing  in  the  x  and 
y  kinetic  energy  expressions  are 
(m3  +  itu  ) 

m  =  — - — - — ,  and 

mA  +  mB  +  mc 


V  = 


IHgltl;. 

mB  +  mc 


which,  for  the  cases  at  hand,  are  listed  in  Tables  2  and  3. 
Notice  that  in  the  limit  where  mA>>(me  +  me)  (e.g.,  for  Ga* 
but  less  so  for  B+)  ,  in  =  me  +  me,  so  bothu  and  in  depend  only 
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on  Che  masses  of  the  B  and  C  atoms,  not  on  the  mass  of  the 
metal  ion. 


II.  b.  Mass  Weighted  Coordinates 

The  above  Hamiltonian,  can  be  rewritten  in  terms  of  mass 
weighted  coordinates  X  =  xV^  and  Y  =  (y  -  yeq)-jjl  as 

„  h2  d2  b2  d2  4  1 

H  =  -Ta?-T5F  +  V(X'Y)  +  5“’Y 

where  V( X ,  Y )  =  v  •  exp(-T(X  -  KY)), 
r  =  — jmw  and 


-  b  m 
a  V/i 


so  TK  = 


are  the  mass  weighted  version  of  the  range  (a)  and  coupling 
(b)  parameters,  and 


cot  =  — 


is  the  square  of  the  Y-mode  harmonic  frequency. 

For  the  ion-molecule  collisions  under  study  here,  the  values 
of  v,  r2,K2  r2#,  and  ci>y  are  given  in  Tables  1,  2,  and  3. 


The  mass  weighted  coordinates  are  introduced  here  and 
used  in  subsequent  classical  and  quantal  calculations  because 
the  findings  of  ref.  3  cause  us  to  emphasize  the  local 
frequencies  of  motion  along  x  and  y  in  anticipation  of  the 


7 


dynamical  resonances  postulated  earlier.  Such  local 
frequencies  are  obtained  by  using  mass-weighted  coordinates 
in  terms  of  which  the  kinetic  energy  is  isotropic  along  all 
coordinates  and  hence  the  potential  energy  contains  all  mass 
and  surface  slope  and  curvature  characteristics. 


II.  c.  Classical  Trajectories  and  Cross-Section 
Evaluation 

We  use  a  classical  trajectory  method  and  our  model 
Hamiltonian  to  compute  the  classical  vibrational  excitation 
cross  section5  as 

=  ]  j  f  f  |[G(y;  v,  J)dy][~dnl[^- d0 

b=C6*G#*Oi)*Oy*y_  L  *  **  JL  ^  7l 

X  [  2  ?rbdb  ]pa(vx,  b ,  v ,  J  ,  y ,  9 ,  <p  ,  rj) 


—  sin  ddd 

JL  2 


where 

Pa  is  the  probability  of  translation  to  vibration  energy 
transfer, 

vx  is  the  asymptotic  relative  speed  of  A  with  respect  to  BC; 
1  -  : 

2  —  ^collision' 

the  angles  <l>  and  ©define  the  BC  initial  orientation, 
r)  defines  BC  rotation  plane, 
b  is  the  impact  parameter, 
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v  and  J  are  the  initial  vibrational  and  rotational  quantum 
numbers  of  BC, 

and  G(Y;V,J)  is  the  distribution  function  for  the  Y 
coordinate  which,  of  course,  depends  on  v  and  J  In  the  cases 
at  hand,  we  take  v=J=0  and  we  approximate  G  in  terms  of  the 
v=0  harmonic  eigenfunction  for  the  unperturbed  Y  coordinate. 

G  ( Y ;  V ,  J )  =  I  <t>0  ( Y )  I2. 

For  our  classical  simulations,  we  focus  on  near-C2v  collisions 
which  therefore  have  small  impact  parameters,  so  we 
approximate  Pa  as 

j,  itrajiv  ,  Y) 

p.  ( Vx,  b,  V,  J,  y,  0,  <f>,  If)  =  5<*i>  )5(tj)4*5(0)5(0  -  -)5(  v)5(J) - - - 

where  ntraj  is  the  total  number  of  trajectories  employed  and 
itraj  (vx,y)  =  1  or  0  depending  whether  the  particular 
collision  has  caused  vibrational  excitation  or  not. 

Of  course,  to  determine  whether  or  not  a  collision  gave 
rise  to  vibrational  excitation,  we  had  to  propagate  the 
corresponding  classical  trajectory.  The  classical  Hamiltonian 
in  mass  weighted  coordinates 

H  =  -P2x+-£^+  -  (ol  Y2  +  v  •  exp  ( -  T  ( X  -  KY ) ) 

2  2  2 

can  be  used  to  obtain  the  following" classical  equations  of 
motion 
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X 


dH 

dPv 


=  P, 


P, 

Y  = 


=  -  — ■  =  V  v  exp ( - T ( X  -  K Y)) 
dX 

dH 


dP, 


=  P, 


dH 


PY  =  -  ~  =  -TK  •  v  •  exp(-nx  -  Ky))  - 
dY 


Q)yY 


We  thus  have  a  set  of  4  coupled  first-order  ordinary 
differential  equations  which  we  chose  to  solve  using  a  fifth- 
order  Runge-Kutta  scheme  monitoring  the  local  truncation 
error  to  ensure  accuracy  and  to  adjust  the  time  stepsize.  We 
used  the  Numerical  Recipes  subroutines  odeint,  rkqs,  rkck4  to 

meet  these  needs. 


For  all  trajectories,  we  set  the  t=0  values  of  the 
coordinates  and  momenta  as  follows 

•  x=10  bohr,  thus  X  =  x-y/ifi  (this  begins  all  trajectories 
where  the  ion-molecule  potential  is  negligible) 

•  Y  and  PY  are  taken  to  be  consistent  with  zero-point 

vibrational  energy  being  in  the  H2  (or  D2  or  HD)  molecule 
11  1 

—  pj  +  —  qjIy2  =  -  h QJV ,  and  Y  is  allowed  to  vary  between  the 
2  2  2 

two  classical  turning  points  with  I  <I>o  (Y)l2  used  as  a 
probability  to  generate  a  series  of  initial  conditions.  For 
each  such  value  of  Y,  both  positive  and  negative  PY  are 

used. 
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the  asymptotic  (incoming)  momentum  along  x  is  determined  by 

Px  ~  —  V2  ^“collision  • 


The  time  duration  of  each  trajectory  is  chosen  to  be 


2 


V  P<  J: 


If  there  were  no  interaction  other  than  a 


specular  reflection  at  X=0,  the  time  for  the  x  coordinate  to 


return  to  x=10  bohr  would  be  -2 


'  x  ' 

V  PX  ) 


.  Because  of  the 


c=0 


repulsive  interaction  potential,  the  actual  time  needed  to 
return  to  x=10  bohr  will  be  less  than  this  estimate.  In 


practice,  we  found  that  by  taking  tf= 


-lU)  , 
2 1  k  L 


x  had 


returned  to  and  passed  10  bohrs  for  all  of  our  trajectories 
and  the  total  energy  was  found  to  be  conserved  within  10“^ 
Hartrees  (i.e.,  to  better  than  0.1%). 


For  each  collision  energy,  a  total  of  200  initial  (Y,Py) 
conditions  were  used.  At  selected  energies,  400  such 
trajectories  were  employed,  but  the  vibration  excitation 
cross-section  changed  little  compared  to  the  200  -  trajectory 
results,  so  the  smaller  number  was  used  in  all  remaining 
cases . 


In  computing  the  cross  section,  itraj(vx,Y)  is  defined  in 

terms  of  the  number  of  trajectories  that  have 
1  1 

-  p2  +  -col Y2  >  h C0V  at  t, ;  this  is  our  criterion  for  defining 
2  *  2  v 

the  Y  mode  to  be  in  v=l  or  higher.  To  probe  excitation  into 
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the  v=l  state  alone,  itraj(vx,Y)  is  taken  to  be  the  number  of 
trajectories  that  have  2h(Uv  >  —  PY  +  —  0)~Y2  >  ho>Y  at  tf. 

Ill.d.  Classical  Trajectory  Results 

Figs.  6a  to  6i  show  the  cross-sections  for  T-V  energy 
transfer  for  MH2+,  MD2+  and  MHD+  for  M=B,  Al,  and  Ga.  These 
classically  evaluated  cross-sections  show  several  noteworthy 
aspects : 

1.  For  B+ ,  Al+ ,  and  Ga+  the  cross  sections  are  indeed 
small  (ca.l0'18cm  )  and  in  the  range  of  what  is  seen 
experimentally.  However,  the  experimental  trend  that  the 
cross-section  is  largest  for  B+  and  smallest  for  Ga+  does  not 
appear  clearly  in  the  classical  results. 

2.  The  T-V  excitation  energy  thresholds  for  B+,  Al+, 
and  Ga'1’  are  "sharp-  and  do  indeed  greatly  exceed  the  reaction 
endothermicities  (which  are  2.6  eV  for  B+  +  H2,  3.9  eV  for  Al+ 
+  H2,  4.1  eV  for  Ga*  +  H2)  and  lie  in  the  neighborhoods  of  3 
eV  for  B+  +  (D2,  HD,  H2),  6  eV  for  Al+  +  (D2,  HD,  H2)  ,  and  9 
eV  for  Ga+  +  (D2,  HD,  H2) ,  not  unlike  the  experimental 
reactive  thresholds.  However,  our  classical  trajectories’ 
prediction  (see  Fig.  6)  that  the  threshold  energies  should 
vary  in  the  order  MD2+  <  MHD+  =  MH2+  is  not  seen 
experimentally  (where  MH2+  =  MD2+*  • 

3.  In  all  M+  +  HD  cases,  there  is,  of  course,  a  single 
classical  T-V  excitation  threshold  energy  as  computed  in  our 
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model.  Experimentally,  MD+  is  formed  at  lower  collision 
energies  than  MH+  when  M*  reacts  with  HD,  so  this  feature  the 
reactive  cross-sections  can  not  be  addressed  in  terms  of  T-V 
energy  transfer  within  our  simple  two-mode  (X,Y)  model.  The 
extension  of  our  analysis  to  treat  the  M+  +  HD  cases  in  terms 
of  three  modes  is  given  in  Sec. III. d. 

III.  Quantum  Treatment 

Although  the  classical  trajectory  numerical  simulation 
data  show  similarities  to  what  is  observed  experimentally,  as 
pointed  out  above,  not  all  features  are  accounted  for.  in 
search  of  reason  underlying  the  remaining  discrepancies  and 
to  gain  further  insight  into  why  the  T-V  thresholds  occur 
where  they  do  and  are  as  sharp  as  they  are,  we  use  analytical 
quantum  tools  on  the  above  model  Hamiltonian.  In  particular, 
we  do  so  to  pursue  a  framework  in  terms  of  which  we  can 
better  explain  the  various  mass  dependencies  of  the  reaction 
cross-sections.  Motivated  by  these  desires  and  again  keeping 
in  mind  the  success  of  the  dynamical  resonance  model 
introduced  in  ref.  3,  we  decided  to  pursue  an  analytically 
soluble,  rather  than  numerical,  refinement  of  the  model 
Hamiltonian  used  thus  far. 

In  ref.  3,  we  made  use  of  local  harmonic  approximations 
to  describe  both  the  H-H  (i.e.,  Y-mode)  and  M+-H2  (i.e.,  X- 
mode)  dependencies  of  the  potential  V(X,Y).  Such  an  approach 
allowed  us  to  specify,  in  terms  of  atomic  masses  and 
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potential  surface  characteristics,  the  dynamical  resonance 
conditions  that  permit  T-V  energy  transfer  and,  by 
assumption,  subsequent  chemical  reaction  to  occur. 

In  this  paper,  we  extend  the  ideas  from  ref.  3  and  use 
locally  defined  harmonic  oscillator  eigenfunctions  as  a  basis 
for  both  the  X  and  Y.  coordinates.  This  may  limit  us  to 
treating  the  early  part  of  the  T-V  energy  transfer  (low  Y- 
mode  excitation)  since,  as  shown  in  Figs.  2  and  5,  the  actual 
potential  energy  surface  is  harmonic  only  in  a  range  of  Y 
values  (ca.  ±  0.4  Bohr)  near  where  T-V  energy  transfer  is 

expected  to  occur.  However,  this  approximation  is  appropriate 
if  the  H2  (D2  or  HD)  is  initially  vibrationally  cold  (as  it 
is  in  the  experiments)  and  if  the  rate  limiting  step  for 
inducing  chemical  reaction  is  the  entrance -channel  T-V 
transfer. 


III. a.  Local  Quadratic  Approximation  to  V 

Focusing  on  the  repulsive  X-dependence  characterizing 
the  entrance-valley  (see  Fig.  2)  part  of  V,  we  define 

F(  X )  =  V( X ,  Y  =  0)  =  v  •  exp(-Tx) , 

and  approximate  F(X)  by  a  local  quadratic  potential,  F(X), 
expanded  in  the  neighborhood  of  a  point  Xo  in  the  region 
where  the  energy  transfer  is  expected  to  occur  from  ref.  3. 
Explicitly, 
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F(  X )  =  Fq  +  A(  X  —  Xn )  + 


1 

'  a‘F 

\ 

2 

X-X0  j 

(X  -  X0)2 


or  equivalently, 

F(  X )  =  F:(l  -  — )  +  -  0)2(X  -  XA)2, 

2  2 

where  the  parameters  in  F  are  chosen  to  make 

32  F  _  d2  F 
<?X2  ”  dx2 


and 

F{X)  =  F( X) 

at  Xo<  thereby  giving 

F0  =  F(X0) , 

and 

=  t2f0  . 

The  choice  of  the  A  parameter,  or  equivalently  of  the 
position  xa  where  the  local  harmonic  approximation  has  its 
minimum. 


where 

Z  =  —  -  — 
<0x  TF0 


implies  a  choice  of  the  derivative  of  F(X)  at  x0. 
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In  independent  work6,  it  has  been  observed  that  equating 
the  fitting  (f)  and  actual  (F)  potentials  at  two  points  often 

produces  a  more  useful  fit  than  is  obtained  by  equating  the 
fitting  and  actual  potentials  and  their  first  derivatives  at 
a  single  point.  Thus,  it  is  not  obvious  that  one  wishes  to 
choose  A  (i.e.,  Z)  to  cause  the  gradients  of  F  and  F  to  be 
equal  at  Xo. 

A.  Fitting  F  to  F  at  Two  Points 

In  terms  of  the  parameters  entering  into  our  local 
quadratic  function  F,  fitting  F  to  F  both  at  Xo  and  at  Xa 
produces  the  following 

F0  =  F{X0), 
and 

F(Xa)  =  F0  ( 1  -  y )  =  F(Xa )  =  ezF0. 

The  second  result  given  above  provides  an  equation 
Z2 

ez  =  (1  - )  that  Z  (or  equivalently  the  slope  parameter 

2 

A=rFoZ)  must  obey.  Solutions  to  this  equation  are  (i)  the 

trivial  solution,  Z=0  (which  does  not  produce  a  second  point 
at  which  F  =  F  since  then  XA=X0)  and  (ii)  Z=-l.  176002,  which 
we  found  numerically.  The  value  Z=-l. 176002, . .is  used 
throughout  the  remainder  of  this  paper  for  reasons  detailed 
shortly. 

B.  Fitting  F  and  F  and  First  Derivatives  of  F  and  F  at  One 
Point 
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If,  alternatively ,  we  require  the  slope  of  F  to  match 
that  of  F  at  the  single  point  Xo,  we  find, 

F.  =  F(X,) 
and 

A  =  -  TF; . 


The  latter  result  implies  Z--1.0  rather  than  Z=-l. 176002. 

A  comparison  among  the  actual  F  potential  and  local 
quadratic  approximations  F  corresponding  to  various  values  of 
Z  is  given  in  Fig.  7a,  and  Fig  7b.  shows  the  difference 
between  F(X)  and  the  various  F(X),  from  which  we  can  clearly 
see  that  z=-l. 176002 .. .  produces  a  potential  that  gives  a  fit 
satisfactory  over  a  larger  interval  of  X.  On  this  basis,  we 
have  chosen  the  local  quadratic  potential  provided  by 


Z=-l. 76002  throughout  this  work. 


In  the  next  Section  of  this  paper,  we  demonstrate  that  a 
resonance  energy  condition  is  met  when  ©x  and  coy  are  related 
by  0)=1.53ci)y.  Thus,  given  (i)y  characteristic  of  H2,  D2,  or  HD, 
we  use  this  identity  to  compute  G)x.  We  then  find  the  point  X: 
at  which  the  model  potential  F=v  exp  (-fX)  has  G)x2  as  its 


second  derivative:  X0  =  -  —  In 


(  p2  \ 

1  V 

\  (0z  , 
\  X  J 


Knowing  Xq  we  can 


evaluate  F(Xo).  The  second  "contact  point"  Xa  is  then  given  as 


X=Xn  -  where  Z  =  -1.176002,  and  the  value  of  F(Xa)  is 

a  u  r 


obtained  as  Fo  ez.  As  also  shown  in  the  next  Section,  the 
point  Xa  is  where  our  model  suggests  efficient  energy 
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transfer  is  to  occur  and  F(Xa)  determines  the  reaction 
threshold. 


This  procedure  then  gives  our  optimal  quadratic  fitting 


function  F 


F(Xa)  +  -0)2(X  -  XA 


in  the  neighborhood  of  Xa-  It  should  be  noted  that  this 
function  matches  the  exponential  F  at  Xo  and  Xa  and  matches 
the  second  derivative  of  F  at  Xo  but  not  at  Xa;  the  second 
derivative  of  F  at  Xa  is  d)2x  while  that  of  F  at  Xa  is  ezQ)x 

(i.e.,  at  Xa  the  local  harmonic  frequency  of  the  actual  F 
function  is  VeF=  0.85  of  its  value  at  Xo.  which  we  express  as 
<o(XA  )  =  0.  85a>*)  . 


Ill.b.  Resonances  in  First-Order  Solution  of  the 
Schrodinger  Equation 


Treating  V(X,Y)-F(X),  which  includes  X-Y  mode  coupling 

and  non-harmonic  X-  character  of  V  as  a  perturbation  Av,  we 

now  express  the  quantum  dynamics  on  the  full  v(x,Y )  +  —Q)2Y2 

2 


2,,2 


surface  in  terms  of  approximate  dynamics  on  the  F(X)+— G>rY 

2 


fit  potential  surface.  We  begin  by  rewriting  the  Schrbdinger 
equation  for  the  2-dimensional  Hamiltonian  as 


iL 

2 


d 2 

dQ 2 


+  F(X.  )  +  -  (i)IQ 2 


—  ~^—r  +  -  Ci)lY2  -  e!'F(Q,  Y)  =  -A  V(Q,  Y)*¥{Q ,  Y 
2  dY2  2  I 


where, 

Q  -  x  -  xA 
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is  the  displacement  along  the  X  coordinate  from  the  point  Xa 
at  which  F= F=Fa,  and  near  where, 

F(Q)  =  Fa  +  |  <02Q2  . 

In  terms  of  the  displacement  coordinates  Q  and  Y,  the 
perturbation  occurring  on  the  right  side  of  the  SchrOdinger 
equation  is  written  as 

-A V(Q,  Y)  =  F(Q)  -  V(X,  Y)  =  FA  +  —  0 )2Q2  -  FAe"riQ_KY) . 

2 

The  eigenvectors  lP(Q,Y)  of  the  full  Hamiltonian 

involving  V(X,Y)  can  be  expanded  in  terms  of  the 
eigenfunctions  { <&„(£>)<!>* (Y)}  of  the  left-hand  side  of  the 

SchrOdinger  equation  as 

n-0  0 


where. 


are  the  conventional  harmonic  oscillator  functions  of  Q  and 
Y,  and  CnM  are  expansion  coefficients. 

To  model  the  entrance -channel  initial  state  of  the  M+-H2 
system,  we  specify  the  quantum  number  (I)  for  the  Y-mode 
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describing  the  H-H  vibration  as  well  as  a  quantum  number  (i) 
relating  to  the  initial  energy  content  of  the  X  mode.  These 
quantum  numbers  thus  specify  an  initial  (unperturbed) 
eigenvector  ( Q ,  Y )  =<1>1  ( Q )  <I>:  ( Y )  ,  and  an  initial  zeroth-order 

energy 

1  1 

.Ej  z  =  Fa  +  ftcox(i  +  — )  +  h(oY(I  +  — ). 

2  2 

Because  we  wish  to  formulate  solutions  'F(Q.Y)  in  terms 
of  a  perturbation  that  induces  transitions  i,  I— m, M  but  not 

energy  shifts,  we  add  to  both  sides  of  the  above  SchrOdinger 
equation  the  average  value  AVu  of  Av  for  the  specific  i,I 

initial  state  of  interest: 

AVn  =  (4>i(Q)0I(Y)  |AV(Q,  Y)|  <&t(Q)Ot(Y)). 


It  can  be  shown  (using  the  expressions  on  p.  60  of  ref.  7) 
that  this  average  value  can  be  expressed  as  follows 


~AVir  =  Fa  +  i(i  +  ±)h(Ox 


FA  e 


^  t 


(k'.)Ui  -  k)'A  2(0 


1 1 

(L  !)2  ( I  -  L)\ 


To  clarify  the  physical  content  of  this  energy  shift, 

r2K2 

expansions  in  powers  of  K2  and  -  can  be  carried  out  (see 

G>y 

table  3  where  it  is  shown  that  these  parameters  are  «1.0)). 
The  lowest  order  terms  thus  obtained  are: 
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~A^jr  =  Fa  +  +  |) 


Recalling  chat  6)*  =  F2F0  and  Chat  FA  /  F0  =  e2  =  0.31  we  find 


So,  the  i,I  diagonal  element  of  Av  produces  a  shift  in  the  X- 

mode  potential  energy  related  to  the  change  in  the  curvature 
along  that  mode  from  (i>\  at  Xg  to  co2(xa)  =  ez(02x  at  XA .  This 

diagonal  element  also  involves  a  (smaller,  because  it  is 
proportional  to  K2)  change  in  the  Y-mode  potential  energy. 

Defining  the  energy  e  relative  to  the  bottom  of  the 
harmonic  potential  (Fa)  plus  the  "shift"  induced  by  AVu 

«  =  E  -  fa-  &vu  , 

we  can  rewrite  the  Schrddinger  equation  with  a  right-hand 
side  that  causes  only  transition  but  no  further  shift  in  the 
initial  state's  energy: 
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TW  +  lalQl  +  =  ( -  AV  +  AV . .  )  *F  (  Q ,  Y ) 


The  Green  Function  of  the  left-hand  side  of  this  SchrOdinger 


equation  is 


G( Q ,  Y ;  Q  1  ,  Y  ’  )  =  £  £  - 


<3>r.(Q)4>*{Y)4>AQ'  )<J>'M(Y'  ) 


n.o  m.o  (n  +  —  )hco„  +  (M  +  —  )h <o,  -  e 
2  2 


and  the  integral  equation  equivalent  to  the  above  SchrOdinger 
equation  and  its  boundary  conditions  is: 

ViQ.Y)  =  4*°  ( Q ,  Y)  +  J  dQ'  dY  '  G(Q ,  Y ;  Q  '  ,  Y  '  )[  AVlZ  -  AV  ](£) 1  ,  Y  ,)4/(  (?  '  ,  Y  '  ). 


Strong  contributions  to  the  above  integral  over  Y'  and  Q*  are 
expected  whenever  E  approaches  +  ^jh(0Y. 


The  relative  importance  of  each  such  "resonance"  (i.e.,  each 

such  n,M  pair)  is  determined  by  the  magnitude  of  the  matrix 
elements  (<&„<!>„  lAv)^^),  and  the  T-V  excitation  threshold  is 

determined  by  the  lowest  value  of  E  at  which  a  strong 
resonance  can  occur. 


We  know  from  our  classical  trajectory  simulations  that 
Av=l  processes  are  dominant  in  the  T-V  excitation,  especially 

near  threshold,  so  we  can  anticipate  that  the  first-order 
correction  to  4*°  obtained  as  the  first  iterate  of  this 

integral  equation 
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¥l(C,  Y)  =  j  dQ'dY  ‘G(Q,  YjQ'  ,  Y'  )[AV„  -  Av](C  *  ,  Y  ')V°  ( Q  '  ,  Y  '  ) 

should  embody  che  primary  effects  for  our  system  in  the 
threshold  energy  regime.  Because  our  zeroth-order 
wavef unction  is  of  the  form 

*F°(Q,  Y)  =  4>1(Q)<D:(Y), 

with  1=0  (because  the  Hydrogen  molecules  are  initially 
vibrationally  cold), the  first-order  wavefunction  correction 
can  be  written  explicitly 
as : 


The  resonance  condition  relating  (Ox  to  %  outlined  in 
the  preceeding  Section  is  obtained  by  requiring  that  Av=l 

single-quantum  T-V  energy  exchange  be  resonant  fcr  the  1,0 
-»0,1  process.  It  can  be  shown  that  the  expectation  value  of 

the  model  Hamiltonian  in  our  basis  is 
Ena  =  {^>n(Q)^M(Y)|H|<I>n(Q)<l>M(Y)) 
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=  (M  +  -)hcoY  +  -  (n  +  -)b(oy  + 
2  2  2  * 


F,ev 


C  M-1  fcT'K*  \ 

,1  4"r  J  , 

\y _ 

n  ! 

f  *ElVl 

M  ! 

,  ( n 

-  j)  !( j  D: 

J 

[l-o  (M 

-  L)  !(L  !); 

m~:k: 

2(0y 


1  11 

Thus,  £,  ;  =  -  b(0Y  +  -  (1  +  -)h(ox 
2  2  2 


+  F  ev  ^  ' 


a  \ 


1  + 


hr 

2(33, 


and 


E0.i  = 


1 

2 


1 

4 


(  xr4  t*rJKJ  \ 

I  tdl.  4<uv  J 

f,  *t2k2) 

i  + - 

l  2(°  {  ; 

The  identity  Eio=Eoi  can  then  be  solved  for  cox  in  terms  of  (Oy, 
with  the  result 


(0 

33 


„  hr2  2 

2  + - e 

4<oy 

l  +  ez 


2 

1  +  ez 


=  1.53 


Thus,  optimal  T-V  energy  transfer  should  occur  at 
geometries  where  <Bx=1.53a>y,  according  to  this  model.  In 


practice,  this  resonance  condition  directs  us  to  seek  regions 
of  the  potential  surface  (Xq)  near  which  G>x=l .  53o>y.  The 


repulsion  energy  Fo  at  this  geometry  then  gives  the  surface's 

range  parameter  T  =  (i)2x  /  F0,  which  allows  the  geometry 
1  176 

XA  =  X0  +  — 1 - ,  where  the  local  potential  surface  curvature 


CO(XA)  =  (0.85)  (1.53)CDy=1.30CDy,  to  be  computed. 


The  relation  between  the  repulsion  energy  FA  at  XA  and 
the  threshold  collision  energy  ECoii  necessary  to  achieve 
facile  T-V  energy  transfer  can  be  seen  by  examining  the 
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energy  dependence  of  che  multitude  of  terms  contributing  to 
vF1(Q,Y).  In  particular,  the  denominator  relating  to  the  0,1 
— »  1,0  process  can  be  expressed  in  terms  of  the  total 


collision  energy  E  as 

1  3 

- h(0  V  +  —  h(l)y  +  F ,  |  1  + 

4  2  * 


ftp 


2  \ 


2(0 


nr2  hK2r2 

e -  + - E , 


X  J 


4  Q)y 


4  (Oy 


Since  a>x  =  1 . 5 


^p2  ^p2j^2 

toy  and  because  - and  -  are  <<  1  (see  Tables  3  and  41 


coY 


C0V 


this  denominator  will  become  small  when  E  =FA.  For  this 
reason,  we  can  predict  threshold  energies  once  F(XA)  is 
obtained. 


For  all  of  the  reactions  studied,  our  primary  findings 
(the  predicted  threshold  energies  FA)  are  summarized  in 
Tables  4  and  5  and  compared  in  Table  5  to  our  classical  T-V 
energy  transfer  onsets  (Eciass),  the  experimentally  observed 
reaction  threshold  energies  (Eexp) »  the  endothermicities 
(Ether)  '  and  the  threshold  predictions  of  our  earlier  work 
( Ecross ) . 

Although  our  FA  thresholds,  which  are  based  on  using  our 
model's  prediction  of  T-V  energy  transfer  onset  as  the 
chemical  reaction  onset,  display  the  general  trends  of  the 
experimental  thresholds  they  still  are  not  in  quantitative 
agreement  and  they  do  not  differentiate  between  MD+  and  MH+ 
thresholds  in  the  HD  case. 
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ZZI.d.  Asymmetry  in  the  MD+/MH+  Thresholds 

0:  e  collisional  energy  has  been  deposited  into  the  B-C 
vibrational  motion,  the  possibility  of  forming  MH+  or  MD+ 
product  ions  is  assumed  to  be  non-negligible  in  our  model. 

For  M++H2  or  M++D2  collisions,  only  MH+  or  MD+ ,  respectively, 
can  be  formed.  However,  in  the  M++HD  case,  a  new  issue 
arises;  it  remains  to  be  explained  why  MD+  is  formed  at 
(significantly)  lower  collision  energies  than  MH+. 

We  begin  our  analysis  of  this  aspect  of  the  reaction  by 
introducing  the  following  three-atom  classical  Hamiltonian 

H=T+U , 

where  the  kinetic  energy  is  written  in  terms  of  the  three 

masses  and  velocities  as 

1  i2  1  1  -2 

1  —  2  mArA  2  mBrB  2 


and  the  internal  potential  energy  is  assumed  to  be  locally 
represented  as  a  sum  of  quadratic  functions  of  the  three 
interatomic  distances 


We  next  transform  to  a  coordinate  system  involving  the 
center  of  mass  coordinate  of  the  three  atoms  ? ^  and  two 

relative  position  vectors  x  and  y  as  shown  in  Fig.  8. 
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Expressing  T  and  U  in  terms  of  these  three  new  vector 
coordinates,  we  obtain 


1  1  i,  1 


u  =  k  ktc  <y  -  ?7c  ?  + 1  (*  -  -7-  y  -  );  +  i  Jc„  (x  +  y  -  r£  f  , 

*  *  M  b  2  M  c 


Where  the  effective  masses  in  and  H  are  as  employed  earlier  in 
this  paper  and  M=mA+mB+mc.  Now  introducing  mass  weighted 
coordinates 

y  =  V£( y  -  Yec.) 

X  =  VA(x  -  x^ ) 


where 


V  =  fe<* 
y*q  '  BC 


and 


?  -  7tq  —  ^  r'q  —  r‘q  4-  ^  rtq 
-Ktq  ~  '  AB  _  'BC  —  'AC  T  _  '  BC 


m 


B 


mr 


We  are  able  to  write  H  in  a  manner  reminiscent  of  the 
Hamiltonian  used  earlier  in  Sec. II,  except  that  x  and  y  are 

vectors  in  a  plane. 

T  =  -X2+  —  Yz+  -  M?2m 
2  2  2  CM 

U  =  -  G72Y2  +  -  GJ2X2  -  SX  •  Y 
2  2 


Clearly,  the  center  of  mass  coordinate  propagates 
independently,  so  it  can  be  removed  from  further 
consideration.  The  5Jy  and  57x  frequencies,  which  need  not  be 
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identical  to  the  G)x  and  0)y  used  earlier  in  this  paper,  are 
expressed  in  terms  of  the  force  constant  parameters  defining 
U  and  the  pertinent  masses  as  follows: 


and  the  parameter  8  governing  the  strength  of  the 
X  and  Y  coupling  is  given  by 

s  =  ^ 

V  m  y  mB  mc 

In  the  m++H2  and  m++D2  cases,  the  8  parameter  vanishes. 


We  can  rewrite  the  potential  as  a  quadratic  form 


U(X,  Y) 


involving  a  matrix 


whose  eigenvalues  are 

(ol  =—  (tfx  +  Oj)  ?  s,  where 
2 
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The  corresponding  two  eigenvectors  of  the  above  matrix  allow 
the  normal  mode  displacement  vectors  of  the  coupled  system  to 
be  written  as 


1  at 

-f  _  4-  _  I  _ I 


2  1®.- 


21®;- 


♦  i-i 


<2  2^®;- 

+  JTTit^I 

1  2 1®;- 


The  AB  and  AC  interatomic  distances  can  be  written  in 
terms  of  X  and  Y  as 


H  1 


AB  _ 


">  f  where 

- =  J 

mc  J 


^AB  ^AB  ^AB 


r  -  req 

AC  •‘■AC 


Then  using  the  above  expressions  for  the  eigenmodes,  these 
same  interatomic  distances  can  be  rewritten  as 


>c  # 


In  terms  of  the  vectors,  the  Hamiltonian  becomes 


h  =  t  +  u  =  -  fit  +  -  a£nt  +  -  fjt  +  -  coiril 


In  the  remaining  portion  of  this  analysis,  we  proceed  as 
follows : 

(i)  We  assume  that  once  the  total  collision  energy  E 
reaches  the  range  F(Xa)  where  appreciable  T-V  energy  transfer 
begins,  the  subsequent  dynamical  evolution  of  the  three-atom 
system  is  best  represented  in  terms  of  time  evolution  of  the 
local  normal  modes 

(ii)  We  partition  the  total  energy  E  in  excess  of  f(Xa) 
between  the  two  normal  modes  E  =  F(Xa)  +  E+  +  E-,  with 

E+  =  rjl  +  0)1  fjl^j  and  E_  =  j  Q)l  describing  the 

energy  content  of  the  two  modes. 

(iii)  We  then  use  equipartition  of  the  excess  energy 
(i.e.,  assuming  adequate  time  to  permit  appreciable 
randomization  of  the  energy  E-F(Xa)  in  excess  of  the  lowest 

threshold)  to  relate  the  mean  square  A-B  and  A-C  interatomic 
displacements  (R2ab)  and  (r\c),  to  (t)*)  and  hence  to  e-F(Xa). 

(iv)  We  show  that,  in  the  BC=HD  case,  (r\h  )  is  greater 
than  (rIq')  for  any  E  >  F(Xa).  This  observation  is  used  to 

infer  that  a  lower  collision  energy  is  required  to  eject  the 
H  atom  (thus  leaving  AD+)  from  this  energized  three-atom 
complex.  Hence  the  threshold  for  MD+  formation  should  occur 
at  lower  total  energies  than  that  for  MH+  formations. 

1.  H2  and  D2  Cases 
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Let  us  first  consider  the  M++H2  and  M++D2  cases  where 


mB=mc=n\H  for  H2  and  mB=mc=2rr\H  for  D2»  and  for  which  8  =  0 .  so 
we  achieve  considerable  simplification  in  the  above 
expressions.  In  particular,  now 


and 

1  _  1 
Vin  2^2n 
1  + _ 1__ 

Va  2J2 ji 

The  average  values  of  (r2u)  and  (r*c)  can  be  expressed  in 
terms  of  those  of  (f}±) 

(rL)  =  <n.  ♦ 

m  \2m  [i 

m  y  2m  n 

and  (r)_  •  V+)  -  0  • 

From  equipartition  of  energy 

=  \s. 

=  ^ho>X“  ■ 
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where  T)*ct  are  quantum  numbers  describing  the  energy  and 

momentum  content  of  the  two  local  normal  modes, 

(specifically,  the  energy  content  above  the  "bottom"  f(Xa)  of 
the  local  harmonic  potential)  we  thus  obtain 

(*L)  =  {*ac)  =  4  hoi-T\-:  + 

m  o  u 


Since  kAB=kAC=k  we  can  rewrite  this  results  as 


{*»)  =  = 


8^j(k+  kac)n 


For  H2  and  D2,  H  =  and  )i  =  mn,  respectively,  and  fh  = 
mA(mB  +  mc) 


mA  +  mg  ■+•  wc 


~  2mH  and  m  =  4  mH , 


respectively.  Hence  for  H2 


V2 hi)lt{ 

V4  ^  & -\j(k  +  kBC )  ' 


while  for  D2 


V8Jan*  8^/(Jc  +  ' 


Thus  for  any  given  energy  E  in  excess  F(XA)  characterized  by 
iliff ,  (k2ah)  will  be  larger  than  (r2m)  by  ca.  41%  (i.e.,  their 

ratio  is  V2  . 


Recall  that  our  classical  trajectory  simulations  showed 
a  somewhat  lower  T-V  energy  transfer  threshold  for  M+  +  D2 
than  for  M+  +  H2,  whereas  the  experimental  reaction  thresholds 
seem  to  be  much  closer  for  D2  and  H2.  The  results  of  the 
preceeding  paragraph  imply  that,  although  M+  +  D2  collisions 
may  lead  to  T-V  excitation  at  lower  E,  it  will  require  more 
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excess  energy  (i.e.,  higher  rj*ee )  to  cause  (r2md)  to  exceed  a 

“critical  bond  breaking"  distance.  For  M+  +  H2,  T-V  excitation 
requires  more  collision  energy,  but  for  to  then  exceed 

the  critical  value,  less  excess  energy  (i.e.,  smaller  )  is 

needed.  The  net  effect  is  that  these  two  competing  tendencies 
essentially  cancel,  thus  rendering  the  H2  and  D2  reactive 
thresholds  very  similar. 

2 .  The  HD  Case 


We  define  mc  =  mH  and  mB  =  mD  =  2mH,  as  a  result  of 

.  .  .  2  kun  5  k  2  2  k  -  -1  l~ in  , 

which  fJ  =  -  mH ,  fi7y  =  ,  +  7  — '  Bx  =  — ,  d  =  — ,  - - Bx 

a  U  n  m..  m  ?.  \  0 


rx  and 


noticing  that  fiJx  <<  GJy  and  S  «  CJy,  we  obtain 


0) 


2 

♦ 


33 


!_  ®Li£ 


2  A 


o>  :  -  ) 


*»C  1M 


"©?-©* 
,  <  -  co: 


(*> 


s2 

Neglecting  terms  in  -3 - r,  oil  *  ©x  and  0)1 


©‘  -  ©x 


©^ .  Thus 


c  1  m 


(  V 

a 

m  g 
K  C  J 


(nl) 


ho)  n 

=  ■■  + 

mco_ 


(  ^ 


V  c  J 


puol 


eft 


More  explicitly. 


■^6kmH 

hT)la 
■^6  kmH 


+ 


+ 


'■BC 


) 


mH 


flT] 


eff 
± _ 


-BC 


+ 


Thus,  for  any  given  excess  energy  E  -  F(Xa)  ( i .  e . ,  ( rje/f ) )  , 
(rah)  will  exceed  (#^). 


The  fact  that  (r*h)  >  (rad)  implies  that  the  AH  bond  is 

stretched  more  than  the  AD  bond  at  any  excess  energy. 
Alternatively,  one  can  conclude  that  less  total  energy  is 
needed  to  "break"  (i.e.,  to  effect  T-V  transfer  and  extend  to 
or  beyond  some  critical  distance)  the  AH  bond,  so  the 
threshold  for  which  H  departs  leaving  AD+  should  lower  than 
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that  when  D  departs  leaving  MH+ .  Indeed,  experimentally,  AD+ 
is  observed  to  form  at  lower  collision  energies  than  AH+ . 

IV .  Summary 

We  have  introduced  a  model  dynamics  to  use  in  treating 
the  T-V  energy  transfer  process  that  seems  to  be  the  rate- 
limiting  step  in  M++  H2  -*  MH+  +  H  reaction,  with  M=B,  Al,  Ga. 

In  applying  this  model,  we 

(1)  fit  our  fully  ab  initio  M+  +  H2  potential  energy 
surfaces  (with  M=B,  Al,  Ga)  to  a  two  dimensional  model 
potential  form,  and 

(2)  Show  how  to  extract  from  local  surface  curvature 
information  the  range  (a)  and  strength  (v)  parameters  needed 
to  use  the  model  put  forth  here  in  a  predictive  manner. 

(3)  Used  this  model  potential  within  a  purely  classical 
trajectory  study  to  conclude  that  collisional-to-vibrational 
energy  transfer  thresholds  seem  to  correlate  reasonably  well 
with  experimental  reaction  thresholds'  for  the  M+  +  H2,  D2,  HD 
cases  at  hand. 

However,  this  classical  treatment  of  the  model  displayed 
significantly  lower  thresholds  for  M+  +  D2  excitation  than  for 
M+  +  H2,  and  it  was  not  of  adequate  detail  to  treat  the  M+  + 

HD  -»  MH+,  MD+  thresholds  asymmetry Therefore,  we 
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(4)  introduced  a  locally  quadratic  approximation  to  the 
potential  surface  to  effect  a  quantal  analysis  of  the  T-V 
energy  transfer  process,  which  suggests  facile  transfer 
occurs  in  regions  of  the  potential  surface  where  certain 
resonance  conditions  are  met. 


(5)  Introduced  a  coupled  three-atom  classical  dynamics 
model  to  examine  internal  mean  square  displacements  in  the 
energized  M+HD  transient  specie-  formed  via  T-V  energy 
transfer,  that  allowed  us  to  suggest  (i)  why  MD+  is  formed  at 
considerably  lower  collision  energies  than  MH+  in  the  M+  +  HD 
reactions,  and  (ii)  why  the  M+  +  H2  and  M+  +  D2  reaction 
thresholds  are  very  similar  although  the  M+  D2  T-V  excitation 
thresholds  are  lower  than  those  for  M+H2. 


(6)  Showed  that  the  probability  of  T-V  Transfer  varies 


as 


nr-K: 


r 


(0,. 


fib2 


,  which  is  a  small  number,  in  agreement 


with  the  small  cross-sections  seen  experimentally. 


In  future  applications,  we  foresee  our  model  dynamics 
being  used  in  either  of  two  modes: 

1.  From  locally  computed  potential  surface  information 
(which  is  assumed  to  be  repulsive  along  one  coordinate  x,  and 
reasonably  harmonic  along  another  y) ,  strength  v  and  range  (a 
or  D  parameters  are  extracted.  Threshold  energies  can  then 
be  predicted  F(XA)  =  ezF0  =  e2  C0x2/r2  =  0 . 122a*?2 m /a2  in  terms 
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of  the  frequency  of  the  BC  mode  to  be  excited,  the  mass  m, 
and  the  repulsive  range  parameter  a. 

2.  Alternatively,  given  experimental  knowledge  of  (T-V 
rate  limited)  reaction  threshold  energies  F  (  Xa  ) ,  one  can 
estimate  the  repuisive  range  parameter  a=  (ez  ( 1 . 53  )  2Ci>y2 
n7/F(XA))1/2  for  various  isotopic  B-C  species  (for  which  m, 
0>y,  and  F(Xa)  vary).  The  same  range  parameter  a  should  be 
determined  for  all  isotopes. 
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Figure  Captions 


Fig.  1  Cross-sections  for  reactions  for  A+=B+(1S)(a)  and 
A+=Al+ ( 1S) (b)  as  a  function  of  kinetic  energy  in  the 
center-of -mass  frame  (lower  scale)  and  laboratory  frame 
(upper  scale) .  Arrows  indicate  the  thermodynamic 
thresholds  for  formation  of  the  X22>  and  A2I1  states  of 

the  products1. 

Fig.  2  (a)  C2v  symmetry  contour  plot  of  the  (xAi)  ground 

state  energy  of  B*  +  H2.  The  R  (the  distance  of  B+  to  the 
center  of  H-H)  and  r  (H-H  distance)  axes  are  in 
Angstroms,  and  the  contours  are  spaced  by  10.0  kcal/mol. 

(b)  C2v  symmetry  contour  plot  of  the  (xAi)  ground  state 
energy  of  Al+  +  H2.  The  R  (the  distance  of  Al+  to  the 
center  of  H-H)  and  r(H-H  distance)  axes  are  in 
Angstroms,  and  the  contours  are  spaced  by  10.6  kcal/mol. 

(c)  C2v  symetry  contour  plot  of  the  (1Ai)  ground  state 
energy  of  Ga+  +  H2.  The  R  (distance  of  Ga+  to  the  center 
of  H-H)  and  r (H-H  distance)  axes  are  in  Angstroms,  and 
the  contours  are  spaced  by  10.4  kcal/mol.  In  (a)-(c), 
the  symbol  X  is  used  to  denote  the  location  of  the 
barrier,  and  Y  is  used  to  denote  the  region  of  strong 
mode  mixing2. 

Fig.  3.  Coordinates  for  near  C2v  -collisions  (A  denotes  the 
M+  ion  and  BC  the  H2,  D2,  or  D)  . 
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Fig.  4. a.  BH2*  Surface  Fit,  where  y=l. 405011  a.u.=.7435 
Angstroms  and  1.75  a.u.  <  x  <  3.5  a.u.  E=39.173 

exp (-2. 81145  x) ;  R2  =  0.936. 

Fig.  4.b.  AlH2+  Surface  Fit,  where  y=l. 405011  a.u. =.7435 

Angstroms  and  2.0  a.u.  <  x  <  3.5  a.u.  y  =  27.271  *  exp 

(-2.091806X) ; R2  =  0.993. 

Fig.  4.c.  GaH2*  Surface  Fit,  rH2=1.405011  a.u.  =  .7435 

Angstroms  and  1.70  a.u.  <  x  <  4.7  a.u.  E  =  14.205  *  exp 
(-1.72899  x) ;  R2  =  0.999. 

Fig.  4.d.  BHD*  Surface  Fit  where  y  =  1.405011  a.u.  =  .7435 
Angstroms  and  1.75  a.u.  <  x  <  3.5  a.u.  E  =  37.182  exp 

(-2.81399  x) ;  R2  =  0.933. 

Fig.  5. a.  BH2+,  x  =  2.50  Bohr,  y*  =  1.460773  Bohr,  k  = 

.369333  a.u.,  E-Evib  =  v  exp  (-ax-b(y-y*) ) ;  In  (E-EVib) 
-2.5809  -  0 . 31132y ;  R2  =  0.847 

Fig.  5.b.  AlH2+,  x  =  2.75  Bohr,  y*  =  1.404289  Bohr,  k  = 

.369333  a.u.,  E-Evib  =  v  exp  ( -ax-b (y-y* ) ) ;  In (E-Evib)  = 
-2.1422  -  0 . 17095y ;  R2=  0.861. 

Fig.  5.c.  GaH2+,  x=3.00  Bohr,  y*=l. 414951  Bohr,  k  =  .369333 
a.u.,  E-Evib  =  exp (-ax-b* (y-y*) ) ;  ln(E-EVib)  =  -2.2229  - 
0 . 18944y ;  R2  =  0.859. 
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Fig.  6. a.  BH2+  Cross  Section.  The  threshold  energy  is  2.8 
eV.  The  grid  is  .5  eV  or  better.  v=39.173  exp(-2. 81145 
x- . 31132  (y-yeq)  )  • 

Fig.  6.b.  BHD+  Cross  Section.  The  threshold  energy  is  2.4  eV. 
The  energy  grid  is  . 5  eV  or  better,  v  =  39.173  exp 

(-2.81145X  -  .31132  (y-yeq) )  . 

Fig.  6.c.  BD2+  Cross  Section.  The  energy  threshold  is  2.4  eV. 
The  energy  grid  is  .5  eV  or  better,  v  =  39.173  exp 

(-2.81145  x  -  .31132  (y-yeq)  )  . 

Fig.  6.d.  AlH2+  Cross  Section.  The  energy  threshold  is  5.1 
eV.  The  energy  grid  is  .5  eV  or  better,  v  =  27.271  exp 

(-2.09181  x  -  .17095  (y-yeq) )  . 

Fig.  6.e.  AlHD+  Cross  Section.  The  energy  threshold  is  5.0 
eV.  The  energy  grid  is  . 5  eV  or  better,  v  =  27.271  exp 

(-2.09181  x  -  .17095  (y-yeq) )  . 

Fig.  6.f.  AlD2+  Cross  Section.  The  energy  threshold  is  4.4 
eV.  The  energy  grid  is  . 5  eV  or  better,  v  =  27.271  exp 

(-2.09181  x  -  .17095  (y-yeq> )  . 

Fig.  6.g.  GaH2+  Cross  Section.  The  energy  threshold  is  5.6 

eV.  The  energy  grid  is  0.5  eV  or  better,  v  =  14.205  exp 

(-1. 72899x-. 18944  (y-yeq)  )  . 
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Fig.  6.h.  GaHD+  Cross  Section.  The  energy  threshold  is  5.7 
eV.  The  energy  grid  is  . 5  eV  or  better.  v=14.205  exp 

( -1.72899x-. 18944 (y-yeq) ) . 

Fig.  6.i.  GaD2+  Cross  Section.  The  energy  threshold  is  5.0 
eV.  The  energy  grid  is  .5  eV  or  better,  v  =  14.205  exp 

(-1.72899X-. 18944  (y-yeq)  )  . 

Fig.  7. a.  F  and  F  for  various  Z,  where  F  =  9.3751  *  exp 
(-1.17600  X),  Xo=l,  XA=2,  Xd=1.8503  and  (Dx-2 . 

Fig.  7.b.  Difference  Between  F  and  F  where,  F(X)  =  9.3751* 
exp  (-1.1760  X)  X0=l,  XD=1 . 8503  ,  XA=2  and  cox=2  . 

Fig.  8.  Internal  Coordinate  System 
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0 


E(Hartree) 


E(Hartree) 


E(Hartree) 


In(E-Euib) 


In(E-Euib) 


0->l  Cross  Section 
Total  Cross  Section 


< 

E 

u 

SO 

l 


0.040 


0.030 


0.020 


U  0.010 

at 

</> 

«i 

2  o.ooo 


■& 

♦ 


5  i  n 

Energy(ell) 


0.050 

0.040 

0.030 

0.020 

0.010 

0.000 


Cross  Section  (10^-16  cnT2) 


-i  -  -  —  >  i  ■  ■  >  >  i — i — ' — ■ — « — i 

0  5  10  15 

Energy(eU) 


0->1  Cross  Section 
Total  Cross  Section 


Cross  Section  (I0A-I6  cm"2) 


0.04 


0.03 


0.02 


0.01 


0.00 


0  5  10  15 

Energy(eU) 


Cross  Section  (10^-16  cm~2) 


Energy(eU) 


nru. 


0-5  10  15 

Energy(ell) 


Function 


Z— 1.176002 
Z—  1 


B+ 

39.173 

exp ( -2 .81145 

x  -0.31132 

(y-yeq) ) 

Al  + 

27.271 

exp (-2.09181 

x  -0.17095 

(y-yeq) ) 

Ga+ 

14.205 

exp ( -1 .72899 

x  -0.18944 

(y-yea) > 

Table  1  .The  v  exp(-ax  -bty-y*,)  )  portion  of  the  U(x,y) 
potentials  used  for  different  M+,  v  is  in  Hartree,  a  and  b 
are  in  Bohr~l. 
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* 


System 

fi  (a .u . ) 

fi)Y  ( a .  u . ) 

H2 

918 

.0200  (.54  eV) 

HD 

1224 

.0174  (.47  eV) 

22 _ 

1836 

.0142  (.39  eV) 

Table  2.  Values  for  the  y-mode  frequency  and  reduced  mass  fd 
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Syst em 

m (a.u. ) 

T (a.u. ) 

1^2  _  m  b 
li  a 

hr:K2  , 

- (10-3) 

CO 

Y 

BH2  + 

3107 

.0504 

.0415 

5.3 

BHD+ 

4328 

.0427 

.0434 

4.5 

BD2  + 

5386 

.0383 

.0360 

3.7 

A1H2  + 

3419 

.0358 

.0249 

1.6 

AlHD+ 

4958 

.0297 

.0270 

1.4 

A1D2  + 

6397 

.0261 

.0233 

1.1 

GaH2  + 

3570 

.0289 

.0467 

1.9 

GaHD+ 

5282 

.0238 

.0518 

1.7 

GaD2  + 

6947 

.0207 

.0454 

1.4 

Table  3.  Collisions!  mass  in,  and  mass-weighted  repulsion  (T) 
and  coupling  (K)  parameters. 
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2 

F0  (a.u. ) 

Fa (a.u. ) 

x0  (bohr ) 

xA (bohr) 

“  i  +  e2 

“A 

CD 

33 

to 

+ 

.0306 
(.83  eV) 

.3686 

(10.03  eV) 

.1137 
(3.09  eV) 

1.6596 
(.878  A) 

2.0779 

(1.099 

A) 

BHD-*- 

.0266 
(.72  eV) 

.3881 

(10.55  eV) 

.1197 
(3.25  eV) 

1.6413 
(.868  A) 

2.0596 

(1.090 

A) 

BD2  + 

.0217 
(.59  eV) 

.3210 
(8.73  eV) 

.0990 
(2.69  eV) 

1.7088 
(.904  A) 

2.1271 

(1.126 

A) 

aih2+ 

.0306 
(.83  eV) 

.7306 

(19.88  eV) 

.2254 
(6.13  eV) 

1.730 
(.916  A) 

2.293 

(1.213 

A) 

A1HD+ 

.0266 
(.72  eV) 

.8021 

(21.83  eV) 

.2475 
(6.73  eV) 

1.686 
(.892  A) 

2.248 

(1.189 

A) 

AlD2  + 

.0217 
(.59  eV) 

.6912 

.  (18.81  eV) 

.2133 
(5.80  eV) 

1.757 
(.930  A) 

2.319 

(1.227 

A) 

GaH2  + 

.0306 
(.83  eV) 

1.1211 
(30.51  eV) 

.3459 
(9.41  eV) 

1.4686 
(.777  A) 

2.1488 

(1.137 

A) 

GaHD+ 

.0266 
(.72  eV) 

1.2491 
(33.99  eV) 

.3854 

(10.48  eV) 

1.4061 
(.744  A) 

2.0862 

(1.104 

A) 

GaD2  + 

.0217 
(.59  eV) 

1.0989 
(29.90  eV) 

.3390 
(9.22  eV) 

1.480 
(.783  A) 

2.1604 

(1.143 

A) 

Table  4.  Values  for  the  x-mode  frequency,  F0,  FA,  x0,  and  xA 
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Species 

fa 

xA 

Fcm 

R  3 

Across 

E  3 

‘-'cross 

AEcher3 

B+  +HH 

3.09 

1.10 

>2.8 

>1.05 

<3.9 

2.6 

3 . 3  ±0 . 1 

B+  +  HD 

3.25 

1.10 

>2.4 

>1 . 00 ; MD 

<4 . 6 

4 . 0±0 . 2 

>1.05; MH 

<3.9 

3 . 0  ±0 . 2 

B++DD 

2.69 

1.13 

>2.4 

>1.05 

<3.9 

3  .  3  ±0 . 1 

Al++HH 

6.13 

1.21 

>5.1 

>1.22 

<6.4 

6 . 6  ±0 . 2 

A1  +  +HD 

6.73 

1.19 

>5.0 

>1 . 16 ; MD 

<7.7 

3.9 

6 . 7  ±0 . 1 

>1.22; MH 

<6.4 

4 . 7  ±0 . 1 

Al++DD 

5.80 

1.23 

>4.4 

>1.22 

<6.4 

6 . 6 ±0 . 1 

Ga++HH 

9.41 

1.14 

>5.6 

>1.21 

<7.4 

N.A. 

Ga++HD 

10.48 

1.10 

>5.7 

>1.15; MD 

<9.0 

4.1 

N.  A. 

>1 . 25 ;MH 

<6.3 

Ga++DD 

9.22 

1.14 

>5.0 

>1.21 

<7.4 

8 . 5  ±0 . 5 

Table  5.  Comparison  between  our  Quantum  and  Classical  model 
predictions  and  those  from  the  mass-weighted  Hessian 
eigenvalues  and  experimental  thresholds .All  the  energies  are 
in  eV,  and  all  distances  are  in  k. 
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